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psj ' Abstract 

In this survey we present the interpretation of isomondromy preserving equations on Rie- 
^ ■ mann surfaces with marked points as reduced Hamiltonian systems. The upstairs space is 

On I the space of smooth connections of GL(N) bundles with simple poles in the marked points. 

We discuss relations of these equations with the Whitham quantization of the Hitchin sys- 
tems and with the classical limit of the Knizhnik-Zamolodchikov-Bernard equations. The 
■ main example is the one-parameter family of Painleve VI equation and its multicomponent 

^\ ' generalization. 
ON ■ 

i ■ 

^ ! 1 Introduction 

The famous Painleve VI equation depends on four free parameters (PVIa^p^-y^s) and has the form 
: d^X_l/l 1 1 \/dXy n 1 1 \ dX 



dt^ 2\X^X-1^ X-t)\dt) \t^t-l^X-t)dt^ 
X{X-l){X-t)f ^t t-1 t{t-l)\ 

It was discovered by B.Gambier [Q| in 1906. He accomplished the Painleve classification program 
of the second order differential equations whose solutions have not movable critical points. This 
equation and its degenerations PV — PI have a lot of applications in classical and quantum 
integrable systems (see, for example 0), topological field theories |^], general relativity ^, 
and in the Seiberg-Witten theory Q. In this paper we discuss two important and interrelated 
aspects of PVI: 

• PVI and isomonodromic deformations of linear differential equations; 

• The Hamiltonian structure of PVI. 

The derivation of PVI equation as the preserving monodromy condition was given by R.Fuchs 
while the Hamiltonian structure of PVI was introduced by J.Malmquist M. 



^Contribution in the Proceedings "Integrability: the Seiberg-Witten and Whitham Equations", Edinburgh, 
September 1998 
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We incorporate the one-parameter family PVI;^2 „2 yi i ^2 =PVIi, in a wide class of non- 

linear equations. They preserve monodromies of systems of linear equations on Riemann curves 
with marked points when the complex structures of curves are changed. These systems come 
from the flatness condition of vector bundles on the curves. We restrict ourself by consider- 
ations of smooth connections with simple poles only, and therefore don't include the Stokes 
phenomena. In such general form the isomondromy preserving equations were considered in . 
Our investigation of these systems is inspired by methods developed in classical and quantum 
integrable systems. In general all the systems can be derived using three different constructions: 
A)The symplectic reduction procedure from free infinite dimensional theory. This approach is 



very similar to the derivation of the Hitchin integrable systems |1C]; 

B) The Whitham quantization of the Hitchin systems 

C) Classical limit of the Knizhnik-Zamolodchikov-Bernard (KZB) equations p^ . 

We discuss these constructions separately and then demonstrate their application on the multi- 
component generalization of PVI,y. The presentation is based for the most part on a previously 



published paper |14|. First, in Sect. 2 we consider the elliptic form of PVI. In this form the 
relations of PVI with the Hitchin systems and KZB equations become transparent. Then we 
discuss these three approaches to the isomonodromic deformations. In Sect. 6 the multicompo- 
nent generalization of PVI is described. Finally, we discussed some open problems related to 
PVI and its generalizations. 

2 Elliptic form of PVI 

1. Elliptization procedure 

Soon after discovering of PVI (|1 . 1[) by Gambler, Painleve presented it in terms of the Weierstrass 



elliptic functions |15]. This paper was almost forgotten for ninety years and the elliptic form 
was rediscovered recently in | p!6| , p^ . We follow the derivation presented in where the 
Hamiltonian form and symmetrices of PVI in terms of elliptic functions are treated. 
Consider the family of elliptic curves 

= C/(Z + Zt) (2.1) 
where r H = {Imr > 0} Let p(n|r) be the Weierstrass function 



p{u\ = -^+Yl (7 — \ — ^2 ~ 7 7 72 ) ' = 

^-^ \[u + mcoi + nuJ2) [muJi + nuJ2j / wi 



(2.2) 



In the most part of the paper we put wi = 1 and p(u[r) = p{u\l,uj2)- p{u\t) uniformize the 
elliptic curve 

Pu{u\t) = 4(p(u|t) - ei(r))(p(u|r) - e2(r))(p(u|T) - e3(r)), (2.3) 

ei = P |r) , (Tq, . . . , Tg) = (0, 1, r, 1 + r). 
We consider two kind of transformations. First one is the lattice action 

u ^ u -\- m + nr, r — > r. (2-4) 
It leaves p{u\t) and Pu{u\t) invariant. The second is the modular transformation by PSL2(Z) 
/ u ,aT + b\ , ,x9 , , ^ / u ,ar + 6^ 



P 



CT + d CT + d 



= (cr + dfp{u\T), pu f ^TTil^^^) = (cr + dfpu{u\T). (2.5) 
/ \CT + d CT + d J 
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Now consider another family of elliptic curves Et B, Y"^ = X{X-l){X-t) parameterized 
by S = {t G \ (0, 1, oo)}. There exists the morphism {Er} {Et} defined as 

{u,r) ^(X= Pm^^Y = ^^,t= '-^) . (2.6) 

V 62 — 61 62 — 61 62 — 61/ 



Theorem 2.1 In terms of {u,t) PVIa,i3,'y,s takes the form 



<fu 1 ^ T- 

— = duU{u\T), ?7(u|t) = — — 2 + yk)' (2-7) 

(ao, ... ,03) = (a, -13,7, ^ - S). 



The proof of the equivalence of (1.1) and (2.7) is based on the Picard-Fuchs equation on 
elliptic curves. The Picard-Fuchs operator 

acting on the holomorphic differential uj = {dE/Bx)/y yields the exact differential 

\dE/B (x-t)'^ ■ '^^^ Picard-Fuchs equation just means that periods of dE/B^/y are annihilated by 

Lt. Using the Picard-Fuchs operator Fuchs proved that PVI ( |1.1D is equivalent to the following 

equation 

t(l - £ ,,,,./y = (o + „^ + y^l^ + - i)i|^) Y. 

The equivalence of (^]^) and ( |2.8[ ) follows from the following equality 

, , lt{t-l)Y d'^X 1 f 1 1 1 \/dX\^ 

til - t)u I dEiB^iy = + ^ - 2 + r^T + r^rj (^J + 

The proof is straightforward. 

Thus, PVI can be written in the form of the so-called //-equation fl^ 



00 



where the right hand side is a special section of the bundle Ef. It can be fixed by the symmetries 
of the equation. 

Under the morphism ( |2.6| ) the holomorphic differential dE/u^ on E^- is transformed in 
dE/B^/y, ^iid ^ ™ Lt- More exactly, the left hand side of ( |2.9| ) takes the form 

dz. 



(61 - 62) (61 - 63) (62 - 61)2 C?^^ - 

Taking into account that 

^ = ^(62 -6i)-3/2 p„(n,r) 
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we come finally to (^.Tp. 

2. Hamiltonian structure 



The hamiltonian form of (|2.7| ) is defined by the standard symplectic form 

= Sv5u, (2.10) 

and the Hamiltonian ^ 

H='^-U{u\t). (2.11) 

Consider the bundle V over the moduli space ^A = i^/PSL2(Z) with the symplectic fibers 
parameterized by the local coordinates {v,u). It plays role of the extended phase space for the 
non-autonomous hamiltonian system (|2.1C| ),( ^rTT ). The equation of motion ( |2.7| ) can be derived 
from the action TonV 

JJT = v5u - H5t. (2.12) 

The symmetries of the non-autonomous hamiltonian systems are determined by the invari- 
ance of the two-form uj on V 



- 5H5t = 5v5u - SH5t. (2.13) 



UJ = to 

It follows from ( p.2[ ) and ( |2.5[ ) that the symmetry group is the semi-direct product of Z + Zr and 
the group r(2) C PSL2(Z). We consider a simplified version of this action in Sect. 7 in detail. 
3. Calogero-Inozemtsev equation and PVI 

Let us introduce the new parameter k and instead of ( |2.13D consider 

w = wo - -SH6t. (2.14) 

K 

It can be achieved by the rescaling the dynamical variables {v,u) and periods uji,uj2 v — 
K~2j y — J. (^-^ —J. ^ ^2. Then, (2/7) takes the form 

= -duU{u\T). (2.15) 
Put r = To + nt^ and consider the system in the limit k — > 0. We come to the equation 



{dt^Y 



-duU{u\To) (2.16) 



corresponding to the autonomous Hamiltonian system with the time-independent potential 
U{u\to). It is just the rank one elliptic Calogero-Inozemtsev equation (C/a,/3,7,5) 0) 0- The 
potential U{u\tq) was considered first by Darboux |2^. It arises also in the soliton theory 
Thus, we have in this limit 

There is the inverse procedure (the Whitham quantization) that allows to construct approx- 
imations of non-autonomous systems starting from integrable autonomous systems. It will be 
discussed in Sect. 4. 

Inozemtsev considered degenerations of U{u\tq) playing with the coupling constants, the 
periods, and u. In this way he obtained the trigonometric, rational and exponential interactions. 
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Presumably, they describe the degenerations of PVI to PV-PI in terms of degenerations of elhptic 
functions. Here is one of his potentials: 

1 1 

ao 2 1" '^1 2 ^ '^2 exp u + as exp 2u. 

sinh u sinh 2u 

In what follows we consider only the subfamily PVIj^ corresponding to aj = . 



3 Isomonodromic deformations 

Here we describe the monodromy preserving equations as reduced Hamiltonian systems. The 
original phase space is infinite-dimensional and almost all degrees of freedom are killed by the 
symplectic reduction. Our approach differs from p^ , where PI-PVI equations are treated as a 
result of symplectic reduction from finite-dimensional space. 

1. Hamiltonian approach 
Let Eg be a Riemann curve of genus g. Consider the space FBuhy^^g of flat vector bundle Vg, 
where G = GL(A^, C) with smooth connection A. The flatness means that its curvature vanishes 

Fj^ = dA+]^[A,A] = Q. (3.1) 

Let us fix the complex structure on Tig. Then for A = (vl, A) we have locally a consistent system 
of matrix differential equations 

{d + A)^ = 0, 
{d + A)-^ = 0. 

We modify this system in the following way. First, introduce formally a parameter k G R {the 
level) and consider the operator kO instead of d in the first equation. Let ;U be a Beltrami 
differential on Tg (/i G n'--'^''^\Tg)). It means that in local coordinates 

/i = n{z, z)-^® dz. It allows to deform the complex structure on such that the new complex 
coordinates are 

/ -N - - / -N de{z,z) 

w = z-e{z,z), w = z, ij{z,z) = - — rr- 

1 — oe[z, z) 

The holomorphic operator d^r, = d + fid is annihilates the one- form dw likewise d annihilates dz. 
We do not touch the anti-holomorphic operator kO. In the new coordinates ( |3.lD takes the form 

F4 = (a + dfi)A - KdA + [A, A] = 0. (3.2) 

Thus, we come to the system 

{Kd + A)^ = 0, (3.3) 

{B + fid + A)"^ = 0. (3.4) 

Represent the Beltrami differential as /x = J2a=i ^afJ^i, where /uj, . . . , is the basis in the tangent 
space to the moduli space AAg of complex structures on S^, (/ = dimA^g = 3(7 — 3, for g > 1, 
I = 1, for g = 1). In other words, t = (ti, . . . ,ti) are coordinates of the tangent vector to Aig- 

To fix a fundamental solution of ( |3.3| ),( |3^ ), impose the following normalization for some 
reference point (zq,zq) G 

^{zo,zo) = I. 
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Let 7 be a homotopically nontrivial cycle in such that (zq, zq) G 7 and 3^ is the corresponding 
monodromy transformation 

y{j) = yij{zo,zo% = Pexp<f A. 

The set of matrices {3^(7)} generates a representation of the fundamental group TTi{Y!,g,zo) in 
GL(A^, C). Independence the monodromy y on the deformations of the complex structure means 
that the linear equations 

day = 0, (a = l,...,0 {da = dtj. (3.5) 



are consistent with (|3.3D , (|3.4|) . 

Proposition 3.1 Equations ( \3. 4 j are consistent with ([S.Sj ), iff 

daA = Q, (a = l,...,0, (3.6) 

^aA = -A^J!i, {a = l,...,l). (3.7) 

K 

The proof is straightforward. 

Proposition 3.2 Equations of motion 6^ ), ^^7^ are Hamiltonian. 
Endow the space FBunj^^c with the symplectic form 

= / <6A,SA>, (<,>=tr), (3.8) 

and the set of Hamiltonians 

Ha = l [ <A,A>fi^^\ {a = l,...,l). (3.9) 

Then (^.6|),( ^?7|) are Hamiltonian equations with respect to uj^^^ and Hg. 

Consider the bundle V over the moduli space A4g with FBunY.,G as the fibers. The triple 
{A, A, t) can be considered as the local coordinates of the total space of the bundle. It is useful 
to consider V as the extended phase space p3|. There is a closed two-form on V 



u; = J^^--Y,6Ha5ta. (3.10) 



Though uj is degenerated on V it produces the equations of motion (3^), (3/7), since the form 
uj^^^ is non-degenerated along the fibers. 

The gauge transformations in the deformed complex structure take the form 

f-^Kdf + f'^Af, A^ f'\d + f,d)f + f'^Af. (3.11) 

The form uj is invariant under these transformations, though its constituents uj^^^ and Hg sepa- 
rately are not invariant. 

Introduce a new couple of the connection components A = {A, A'), where A' = A — ^fiA. In 
terms of {A, A') the form to ( p. 10 ) takes the canonical form 



<6A,5A'>. (3.12) 
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2.Symplectic reduction 

The form to is degenerated on FBunj^^Ci because it is invariant under the action of the group 
Q of gauge transformations ( p. Ill ), generating by the flatness condition ( |3.2D . The gauge fixing 
along with the flatness condition (3.2) is nothing else as the symplectic reduction from the space 
of smooth connections Sm^^c in the bundle Vq to the reduced space 

FBun^^C = FBunY,,G/Q = SmY.fil jQ 

The double slashes means that we impose the moment constraints ( p.2[ ) and fix the gauge. 
FBuriY^^G is the moduli space of Eat connections of the bundle Vq- In terms of the symplectic 
reduction procedure the flatness condition is called the moment constraint equation. 

Let us fix the gauge in a such way that the A component of A becomes anti-holomorphic 



dL = 0, {L = f-\d + fid)f + f-'Af). 



(3.13) 



We can do it because the antiholomorphity of f^^{d + fJ-d)f + f^^Af amounts to the classical 
equations of motion for the Wess-Zumino-Witten functional Swzwif, ^) for the gauge field / 
in the external field A. Denote the gauge transformed field A as L 



Then (^]^) takes the form 



L = r^Kdf + r'Af. 



(9 + ^^l)L + [L, L] = 0. 



(3.14) 



Thus, the moduli space of flat connections FBunj^^c are characterized by the set of solutions 
of the linear differential equation ( ^.14 ) along with the condition ( ^.13 ). The moduli space 
FBun-^^G is finite-dimensional space 

dim F^ns,G = 2(iV^ - 1)(5 - 1), g > 1. 

After the gauge fixing we come to the bundle V over Mg with FBuhy^^g as the fibers. The 
system of linear differential equations ( |3.3D ,( |3^ and ( |3.5| ) after the gauge fixing takes the form 

(k(9 + L)^' = 0, (3.15) 

(a + /i5 + Z)^ = 0, (3.16) 

(kS, + M,)^ = 0, (3.17) 

where we replaced ^ on f~^^ and Ms = ndsf f~^- 

The gauge transformations do not spoil the consistency of the system. The consistency (3.15) 
and ( |3.16|) is provided by ( ^.14 ) and ( 3.13|) . In fact, the consistency ( 3.17 ) with ( 3.15 ) and ( 3.16 ) 
leads to the Lax form of the equations of isomonodromic deformations 

dsL - KdM + [M, L] = 0, (3.18) 

KdsL - ={d + fid)Ms - [Ms, L]. (3.19) 

They play the role (p.6|),( |3?7| ) correspondingly. The last equation allows to find Ms in terms of 
dynamical variables L,L. 

The symplectic form uj on the reduced phase space V is 

1 



UJ 



I <dL,6L> --Y6Hs6ts. 



(3.20) 
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H, = l[ <6L,6L>fif^ (3.21) 

Introduce the local coordinates (v,u) in FBun-^^G- 

L = L(v,u, t), L = L(v,u, t), 
V = (7;i,...,t;(;v2-i)(9-i)), u= (Mi,...,M(^2_i)(g_i)). 
Assume for simplicity that this parameterization leads to the canonical form on FBunj^^c 

u;(0) = [ < 6L{v, u, t), 5Z(v, u, t) >= ((^v, 6vi), (3.22) 

where the pairing in the right hand side is induced by the trace. The form on the extended 
phase space is 

w = {6v, 6u)-- Y^ (v, u, t)6ts, (3.23) 

s 

and the variations of the Hamiltonians Kg in the new variables take the form 

SKs = f [<L,6L> + k(< dL, dsL> - < dsL, 6L >)]. (3.24) 



Now, due to (3.14), the hamiltonians depends explicitly on times. Consider the one-form (the 
integral invariant of Poincare-Cartan) 

^ s 

There exist 3g — 3 = dim A4g-dimensional space of vector fields Vg that annihilates 6 

Vs = Kds + {Hs,-}, {s = l,...,l). (3.25) 

It can be checked that Vs satisfy the following conditions 

KdsHr - KdrHs + {Hs,Hr}^m = 0. (3.26) 

Thereby, they define the flat connection in the bundle V. These conditions are called the 
Whitham hierarchy (WH). The equations for any function /(v,u, t) on V take the form 

4f(v,u,t) df{v,u,t) 

They are called the hierarchy of isomonodromic deformations (HID). 

The both hierarchies can be derived from variations of the prepotential J^. It is defined as 
the integral over the classical trajectories in the extended phase space P 

T{u,t) = ^(uo,to) + Csdts, (3.28) 

Juo,to 

where >Cs(5sU, u,t) = {v,dsu) — Ks{v, u,t), (S^u = ^^) is the Lagrangian. T satisfies the set 
of the Hamilton-Jacobi equations 

ndsT + Hs{—,u,t) = 0. (3.29) 
ou 
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The logarithm of J- is cahed the tau-function of HID. 
3. Singular curves 

The singular curves are important for applications, since they produce nontrivial systems 
for the low genus curves {g = 0,1). In these cases the explicit calculations of hamiltonians are 
available. 

Consider a curve ^g,n of genus g with n marked points (xi, . . . , Xn). The number of times is 
equal to dimension of the moduli space Mg,n- We extend the space of connections FBunY,,G = 
{A, A} by adding the coadoint orbits of G = GL(A^, C) in the marked points 

(01,..., a), 0, = {p, = gplg-^}, 

where fixes the conjugacy class of Oh- We allow the A component of connection to have 
simple poles at the marked points, while the Beltrami differentials vanish there. Then instead 
of (^) we obtain 

n 

{B + d^i)A- KdA + [A, A]=Y,^'^ {xb)Pb (3.30) 

6=1 

The Hamiltonian formalism is provided by the modified symplectic form 

= / <5A,5A> +2^2 < 5{pi,g^^),8gi, > . (3.31) 

^9." 6=1 

To derive HID one should start from the space of connections with simple poles in the marked 
points. 

Finally, we come to the same hnear system (Isl^ ), (|3T^),( ^) , but due to (|33c| ) the 
following relation between L and L holds 

n 

{d + dfi)L+[L,L]=Y,S\xb)pb. (3.32) 

6=1 

As before, the linear equations are equivalent to the equations of motion of HID coming from 
the symplectic form 

1 ' 

= cj(o)(v,u,p) - - V5ir^(v,u,p,t)5ts, {I = dim{Mgn), (3.33) 
where p = {pi, . . . ,p„),a;o is determined by the reduction from ( |3.31| ) 




< 6L{^r, u, p), (5L(v, u, p) > +27ri ^ < 5(pftC/^ 6gb >, (3.34) 

6=1 



and K, ( p^ ). 

4 Hitchin systems and their Whitham deformations 

1. Hitchin systems 

Consider the moduli space TZg^N of stable holomorphic GL(A^, C) vector bundles V over T,g. It 
is a smooth variety of dimension 

dimng,N= 9 = N\g- + (4.1) 
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Let T*TZg^N be the cotangent bundle to TZg^N with the standard symplectic form on it. Hitchin 
|10| defined a completely integrable system on T*TZg^N- 

The space T*TZg,N can be obtained by the symplectic reduction from the space T*TZg ^ = 
{^,A), where ^ is a smooth connection of the stable bundle corresponding to d + A and $ is 
the Higgs field $ G Q^{Tig, EndV (S) K) (K is the canonical bundle of S^). There is the well 
defined symplectic form on this space 



(0) = / < S^,6A > . (4.2) 



This form is invariant with respect to the gauge group Q = C^Map{Tig, GL(A^, C)) action 

$^ /-!<!)/, A^ f-'df + r'Af. (4.3) 

In particular, TZg^N = T^gN/^- Ps,k = Ps,kdz~^ ® dz be the {—k + 1, l)-differentials {ps,k ^ 
H^{Tjg,T^^^ (g) K)), and s enumerates the basis in ff^(Sg, F'^"-^ (g) K). {ps,2 = Ps)- Due to the 
Riemann-Roch theorem 

dimF^(Sg,r^-^ (8)ir) = {2k - l){g - 1). 
These differentials allows to define the gauge invariant Hamiltonians 

Hs,k = T [ <'^' > Ps,k, {k = l,...,N, s = l,...,{2k-l){g-l)). (4.4) 

The Hamiltonian equations take the form 

da<^> = 0, {da = ^,a = is,k)), (4.5) 

Ota 

daA = ^'^-Vs.fe (4.6) 

The gauge action produces the moment map p : T*TZ^ Lie*{GL{N, C)). It follows from 
(|4.2|) , (|4.3| ) that p = + [A,^]. The reduced phase space is the cotangent bundle we started 
with 

T*ng,N ~ T*ny/g ■.= p-\o)/g. 

The Hitchin hierarchy (HH) is the set of Hamiltonian equations with Hs i^. (4.4) on the reduced 
phase space T*TZg^N- Hitchin observed that the number of integrals Hg }^ 

N 

Y^(2k-l){g-l) = N\g-l) + l 
k=l 



coincides with dimension g of the coordinate space TZg,]\f Since they are independent and 

Poisson-commute, HH is the set of completely integrable Hamiltonian systems. 
Let fix the gauge of the field A 

A = fdr' + fLf-\ 

Then 

is a solution of the moment constraint equation 

dL+[L,L]=0. (4.7) 

The space of solutions of this equation is isomorphic to H^{T,g, End V <8> K) - the cotangent 
space to the moduli space Tlg^N- 

The gauge transformation / defines the element Ma in Lie(GL(A^, C)) Ma = daff^^- 
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Proposition 4.1 The system of linear equations 

{X + L)Y = 0, 
{d + Y,X''-\kPs,k + L)Y = 0, 

s,k 

{da + Ma)Y = 

is consistent and defines the equations of motion for HH. 



(4.8) 
(4.9) 

(4.10) 



Proof. The consistency of ( [4.8D and ( [4. Of ) follows from (|4.7] ). In terms of L the equations of 
motion ( [4.5| ),( |4^ ) take the form 

daL + [Ma, L] = (the Lax equation), (4.11) 

daL-dMa + [Ma,L]=L''-^Ps,k, {a = {s,k)). (4.12) 

The Lax equation provides the consistency of ([4.8| ) and ( l.lOf ), while the second equation of 
motion ( [4.12 ) plays the same role for the couple ( [4. 9) ) and (4.10). This equation allows to 
determine Ma from L and L. 

The bundles over singular curves can be incorporated in this approach as well |^]. The 
Higgs field has simple poles at the marked points and (4.7) is replaced by 

n 

dL+[L,L] = 2TTiJ2S\xb)pb- 

6=1 

The form on T*TZg^N 

. n 

.(0) = / <6L,5L> +2T,i < S{pbgb'),S9b > 

■'^9 6=1 



(4.13) 



(4.14) 



Note that the space T*lZg,M and the space of flat bundles FBun-^^c have the same dimen- 
sions. Moreover, it follows from the comparison the moment constraints ( |3.32| ),( p.l3 ) and the 
symplectic forms ( |3.34| ),( ^.14 ) that they are isomorphic as symplectic manifolds. 



2. Spectral description 



Due to the Liouville theorem the phase flows of HH are restricted to the Abelian varieties, 
corresponding to a level set of the Hamiltonians k = Cs^k- The phase flow takes a simple 
form in terms of action-angle coordinates. They are defined in a such way that the angle type 
coordinates are angular coordinates on the Abelian variety, and the hamiltonians depend on the 
action coordinates only. To describe them consider the characteristic polynomial of the matrix 
L 

P(A, z) = det(A + L) = A^ + ^lA^^^ + . . . + bjX^-^ + . . . + (4.15) 
bj = Y^Miuj, (Miuj — principle minors of order j, hi = trL, bjsf = detL). 
The spectral curve C cT*T,g is defined as the zero set of P 

C = {P(A,z) =0}. 

It is a well defined object, because the coefficients bj are gauge invariant. 
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Since L G H^{Ylg, End V (E> K), the coefficients bj E H^{Tig, K^) and we obtain the map 

p : T*lZg^N ^B = (Bf=iH\^g, K^). (4.16) 

The space B can be considered as the moduh space of the family of spectral curves parameterized 
by the Hamiltonians Hg^k- The fibers of p are Lagrangian subvarieties of T*TZg^N- The spectral 
curve C is the A^-fold covering of the basis curve S^^tv 

vr : C — > Sg,Ar. 

Its genus g{C) is equal to dimension g of TZg.N- There is a line bundle C with an eigenspace of 
L(z) corresponding to the eigenvalue A as a fiber over a generic point (A, z) 

C C ker(A + L) C 7r*(y). 

It defines a point of the Jacobian Jac{C), the Liouvillean variety of dimension g = g{C). 

Conversely, if z G Eg is not a branch point one can reconstruct V for a given line bundle on 
C as 

Let tOj, j = 1 . . . ,g be the canonical holomorphic one-differentials on C such that for the 
cycles ai, . . . ,ag] j3i, . . . ,Pg, Oi ■ aj = (3i ■ Pj =0, Oi ■ f3j = 6ij, ivj = 5ij. Then the 
symplectic form (PD can be written in the form 



N 

(0) 



/ < 6L,5L >=Y. f ^^jK 



i=i 



Here are diagonal elements of sLs ^, where s diagonalizes L, sLs ^ = diag(Ai, . . . , Aat). 
Then we obtain 

Because A is a holomorphic one- form on C, it can be decomposed as A = X!f=i (^j^j- Thereby 

.7 = 1 



The action variables can be identify with 

Xujj, (j = l...,5). (4.17) 
To define the angle variables, put locally ^ = dlogip. If {pm) is a divisor of ip then 

JC ^ Jvn 



Thus ifj are linear coordinates on Jac{C) and 
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,(0) = 



w^"' = ^ 6aj6(pj. 
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3. Scaling limit. 

Consider HID in the limit k — > 0. The value k = is called critical. We prove that on the 
critical level HID coincide with part of HH relating to the quadratic Hamiltonians ( |4.4| ) . 

Note first, that in this limit the A-connection is transformed in the Higgs field A — > 
and therefore 



FBunj: G ^ T*n. 



9,N- 



But the form cv on the extended phase V appears to be sing ular (see (P^, (l3!20|) ). To get 
around we rescale the times 



t = T + Kt 



H 



(4.18) 



where are the fast (Hitchin) times and T are the slow times. Assume that only fast times 
are dynamical. It means that 



After this rescaling the forms ( |3.10| ),( p.2O ) become regular. The rescaling procedure means that 



we blow up a vicinity of the fixed point /is ^ in Mg^n and the whole dynamic is developed in this 

(4.19) 



vicinity. This fixed point is defined by the complex coordinates 

wq = z -^Tses{z,z), wo = z. 



Now compare the Baker- Akhiezer function of HID ( p^ , (|3l^ ),(p7| ) with the Baker- 
Akhiezer function of HH Y ( |4.8| ),( |4lg| ),( [4.10D . Using the WKB approximation, assume that 



5(0) 

^ = $exp( 

K 



(4.20) 



where is a group valued function and S^^\ S^'^'^ are diagonal matrices. Let substitute ( [4.20| ) 
in the linear system (l3l^),(|3l^),(|3T7|). If 



there are no terms of order k ^. It follows from the definition of the fixed point in the moduli 

(4.21) 



of complex structures ( 4.19 ) that 

5(o)=5(o)(ri,...,r;|z-^r,e,(z,z)). 



We take also S^^'^ = (95(°) Y.s t?^s 

{z,z). In the quasi-classical limit we put 
9^(0) = A. 



(4.22) 



In the zero order approximation we come to the linear system of HH ( |3.15| ),( ^.16|) , (13.171) , 
defining by the Hamiltonians g, k = 1, 2. The Baker-Akhiezer function Y takes the form 



(4.23) 



Our goal is the inverse problem. We need to reconstruct the dependence on the slow times T 
starting from solutions of HH. Since T is a vector in the tangent space to the moduli of curves 
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■Mg^n, it defines a deformation of the spectral curve in the space B (see (p|), (p^ ). Solutions 



Y of the linear systems (^) ( [4.9[) ,( 4.10| ) take the form Y = ^e^^*" ^% where Qs are diagonal 
matrices. Their entries are primitive functions of meromorphic differentials with singularities 



matching the corresponding poles of L. Then according with ( 4.23 ) we can assume that 



_d_ 



■dS = dQc,. 



These equation define the approximation to the phase of ^ in the linear problems (3. 15), (3. 16) 
iKlf) of HID along with 

d 



dS 



The differential dS plays role of the Seiberg-Witten differential. Important point is that only 
part of the spectral moduli, connected with -fffc^s, k = 1,2, is deformed. As a result there 
is no matching between the action parameters of the spectral curve o^, j = 1, . . . ,^ ( 4.17 ) and 
deformed hamiltonians. The detailed analyses of this situation in the rational case is undertaken 
in PI. 



Another object of the Whitham quantization is the prepotential JF (|3.28| ). It depends on the 
action variables aj . This dependence is compatible with the Hamilton-Jacobi equation ( |3.29| ) 
with slow times Tg as the independent variables. These equations are discussed in [E^, [2^. 



5 Classical limit of the Knizhnik-Zamolodchikov-Bernard equa- 
tions 

The Knizhnik-Zamolodchikov-Bernard equations (KZB) are the system of differential equations 
having the form of the non-stationer Schrodinger equations with the times coming from n 
(see, for example, p7[|). 



They arise in the geometric quantization of the moduli of flat bundles FBun-^^c 




Let V = Vi X ■ ■ ■ ^ Vn he the tensor product of finite-dimensional irreducible representations 
associated with the marked points. The Hilbert space of the quantum system is a space of 
sections of the bundle Sy ^quant (Yj g ^Yi) ovcr F BuTi'^ Q depending on an negative number 
with the y-fibers. It is the space of conformal blocks of the WZW theory on Sg^„. 

The Hitchin systems are the classical limit of the KZB equations on the critical level p^ . 
The classical limit means that one replaces operators by their symbols and generators of finite- 
dimensional representations in the vertex operator acting in the spaces Vj by the corresponding 
elements of coadjoint orbits. To pass to the classical limit in the KZB equations 

(k''"'^"*^, + ^,)F = 0. (5.1) 

we replace the conformal block by its quasi-classical expression 

F = expp (5.2) 

where h = (^ij^q"*)-!^ Consider the classical limit ^9"°"* — > oo and assume that values of the 
Casimirs C*, {i = 1,. . . ,rankG, a = 1, . . . ,n) corresponding to the irreducible representations 
defining the vertex operators also go to infinity. Let all values lim ^^uant are finite. It allows to 
fix the coadjoint orbits in the marked points. In the classical limit ( [5.1D is transformed to the 
Hamilton-Jacobi equation for the action = logr (|3.2S| ) of HID. 
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The integral representations of conformal blocks are known for WZW theories over rational 
and elliptic curves |3^, 34]. Then ( ^I^ ) allows to extract the prepotential T of HID. 

The KZB operators ( ^.1| ) play role of flat connections in the bundle P"?"""* over the moduli 
of curves M.g^n with the fibers f y ,,guant |3^ 

+ H,, K9"'*"*a, + H,\ = 0. 
These equations is the quantum counterpart of the Whitham hierarchy ( 3.26| ). 



6 Multicomponent generalization of PVIj, 

Consider Fi?uns,G over the family of elliptic curves with a one marked point The space 

A^i^i is one-dimensional, because the position of one point on a torus is irrelevant. Thus, we 
have only one time r and Adi^i ~ Ej- (|2.l| ). In this case the Beltrami differential takes the form 

T -To 

fJ- = ^• 

T -To 

Consider the most degenerated orbit O = (gp^g^^) of GL(A^, C) sitting in the marked point 
z = with 

= u[{l^_^f (g, - Id]. (6.1) 

N N 

For stable bundles the gauge transforms allow to put A component in the diagonal form 

27ri 

L = ^u, u = diag(ni, . . . , un) G H — Cartan algebra. (6-2) 

r - To 

It means that 

Ldwdw = u. (6.3) 



Let L = d log (p. Then the integral 

i-P 

lEr JPo 



/ Ldwdw = / log (j)dw. 
Je^ JPn 



defines the Abel map Ej- in the product of Jacobians. 

The remaining gauge transforms do not change the gauge fixing. These transformations are 
generated by the Weyl subgroup of G and elements f{w,w) G Map(r2, Cartan(G)). The 
orbit variables can be gauged away by these transforms and we are left with The 
solution L of the moment constraint 

dii,L + [L,L\ = 27rz52(0)pW 

takes the form 

L = P + X, P = 2Tii{— (6.4) 

1 -/i p 

u = diag(ni, . . .,un), v = diag(vi, . . .,vn) 

_ w — w 

Xjk = x{uj - Uk) = (r - ro)z^exp27rz{ —{uj - Uk)}(t){a{uj - Uk),w), 

T - To 
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o{u)e{z) ^ 

The operator M defining the phase flow according with the Lax equation ( |3.18| ) can be extracted 
from (|3l9| ) 

JV ^ 
M = —D + Y, D = diag((ii, . . . ^dj^), dj = ^ s{uj — Ui), s{u) = —p{u) + const. (6.5) 

. , . K 

Yjk = y{uj - Uk), y{u,w,w) = . / — —dux{u,w,w). 

■' 2TnK{T — To) 

The functions x,y,z satisfy the functional equations 

x{u, z, z)y{v, z, z) — x{v, z, z)y{u, z, z) = {s{v) — s{u))x{u + v, z, z). (6.6) 

This equation is derived from the Lax equation ( ^.ISf ). 
The symplectic form lv is boiled down to 

u) = ((5v, (5u) 6H6t, 

where 

(5v,5v) ^"^ ^ ( IN 
-n = -, TTT > P\Uj — Uk\T). 

^ ' j<k 



They define the hamiltonian flow 



^ ' k<j 

FoiN = 2 one can put ui = —U2 = u. Then the potential 



^Pu{Uj - Uk\T). (6.7) 



(2vri)2 ° ^ ' ' (2vri)2 f-^^ 



produces PVIi/ (see (2.7). 

The remaining gauge symmetries implies that uj is invariant under the Weyl transformations 
W of (v, u) and the lattice actions (compare with (E 



V — > sv, V + Kn, u — > su, u — m + rn, (s G n G Z'^). 

It is also invariant under the PSL2(Z) action on r 

ar + b , , ,,1 

r — > -, V v(cr + a) — kcu, u u(cr + a) . 

cr + a 

This invariance is follows from the invariance of the upstairs system under the diffeomorphisms 
of ^g,n (see iQ for detailes). 

On the critical level k ^ 0, t — tq = nt we obtain the elliptic Calogero N-body system. This 
system is a particular example of the Hitchin systems |24]. Note, that the functions x,y and s 
defining the Lax matricies satisfy the same functional equation ( |6.6D as in the Calogero-Hitchin 
limit K = 
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7 Conclusion 

Here we propose a few open problems in the context of topics discussed above. 



• The evident problem is a description of PVI with four arbitrary constants as a reduced 
Hamiltonian system. The first step in this direction is the Lax form of PV I^^p^^^s- The Lax 
form is even unknown on the critical level, i.e. for the Calogero-Inozemtsev system. It will 
be interesting to generalize this approach to the A^-body Calogero-Inozemtsev system and the 
A^-component PVI with four coupling constants. 

• The degenerations of PVI to PV-PI in terms of elliptic functions. 



• There exists a generalization of the Calogero systems related to any simple group. In 
addition to degrees of freedom coming from the moduli of bundles (the coordinates of particles), 
these systems certainly contain degrees of freedom related to the coadjoint orbits. Recently a new 
Lax equations based on arbitrary root systems without the orbit coordinates were proposed 
37 1 . This construction is purely algebraic and does not use the symplectic reduction. How these 



systems can be incorporated in the Hitchin approach, or, more generally, in the isomonodromic 
deformation construction? 

• Consider the N = 2 elliptic Calogero system. The solution n(t), corresponding to the fixed 
value /i2 of the Hamiltonian 

^ = y + 4:^P(2^iko) = ^2, 
is implicitly described by the elliptic integral of the first kind 



1 /"^u fix ^ ^ 
i2uo y' 



where y = 4(x — ei(ro))(x — e2(ro))(x — e3(ro)). As it was mentioned at the end of Sect. 6 it 
can serve for the calculations of solutions to PVIj^. This procedure can be accomplished by 
the Krichever averaging method |3^. It will be interesting to compare this approximation with 
explicit solutions presented recently in ^ for some particular value of the coupling constant 
u. 

As suggested in Sect. 5 another way of approximation comes from the classical limit of confor- 



mal blocks for SL2{C) theory on elliptic curves with one marked point [33, 4C]. Which method 
gives the better approximation? 



• We considered deformations with respect to the moduli of complex structures of curves. 
They describe only part of the moduli of the spectral curves C. The remaining moduli of C come 
from /9<j fc, k > 2. They correspond to the so-called W-geometry of the basic curve ^g,n- This 
geometry is poorly understood. On the other side, there are no examples of isomonodromic 
deformation equations with respect to these moduli spaces, as well as the corresponding higher 
order KZB equations. Any progress in understanding of one of these subjects will shed light on 
another. 
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